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I. INTRODUCTION 


A meeting of the committee for the Extension to the U.S. Standard 
Atmosphere (COESA) , was held in Boston on September 13, 14, and 15, 1971. 
At this meeting it was agreed to review the current state of knowledge 
of the earth's atmosphere with the objective of possibly revising the 
"U. S. Standard Atmosphere, 1962" above the altitude of 50-60 km, and 
of extending the upper limit of the tabulations of a revised standard 
to 1000 km. 

Various schemes have been used to divide the vertical extent of 
the atmosphere into a number of regions. One of these is based upon 
the major features of the temperature-height profile. Another is based 
upon major features in height distribution of the constituent gases of 
the atmosphere. In order to facilitate the review of data and the re- 
vision of the "U.S. Standard Atmosphere, 1962" C0ESA has chosen to 
divide the atmosphere in accordance with a two-fold criterion: 

(1) the principal techniques used for measuring the thermodynamic 
properties of the atmosphere (temperature, density and pressure) and 

(2) the amount of observational data available for constructing 
and testing any proposed model. Applying this criterion of measuring 
technique and quantity of data within the height region of interest, 

50 to 1000 km, COESA has divided the atmosphere into three overlapping 
regions, and has assigned a particular task group to concentrate on 
each of these three regions : 

Task Group I, 50 to 100 km 

Task Group II, 80 to 200 km 
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Task Group III, 140 to 1000 km 

The height region for Task Group I, 50 to 100 km, is a region in 
which atmospheric measurements of temperature, density and pressure 
are made almost exclusively with rocket -borne instruments; these have 
served to develop an extensive set of thermodynamic data. The height 
region for Task Group III is one in which the thermodynamic properties 
are determined almost exclusively from satellite related observations: 

A vast amount of data, particularly mass-density data has been accumulated 
for this height region, 140 to 1000 km. 

The height region for Task Group II, primarily the 40 kilometers 
between 100 and 140 km, plus an overlap into the regions of the other 
two task groups is a region for which niether the rocket-related observa- 
tions nor the satellite-related observations provide much information. No 
unique vehicle or observational technique has to date been developed for 
efficient observation of this height region, and only a limited amount of 
thermodynamic data for this region were available for consideration by Task 
Group II; consequently, the job of Task Group II became primarily one 
of fitting an analytical function to bridge the gap between models which 
were appropriate to the other two regions. 

Task Group II met on four different occasions 

February 22, 1972 GSFC Greenbelt, Maryland 

April 19, 1972 Sheraton Park Hotel, Washington, D.C. 

November 3, 1972 GSFC Greenbelt, Maryland 

December 8, 1972 AFCRL Bedford, Massachusetts 

with Champion, Reber and Minzner in attendance at each of these meetings 
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and with other members absent from one or more meetings. Only two 
basic proposals were submitted to the committee, one by G. R. Swenson, 
and one by C. A. Reber. The first two meetings were devoted to the 
discussion of these proposals, with successive revisions between meet- 
ings. By the third meeting, some of the differences between these two 
proposals had been resolved and Swenson withdrew his proposal in favor 
of the successively revised Reber model. At this meeting empirical 
number-density concentrations of various atmospheric species were agreed 
upon by the committee for the 150-km height level. Some revisions to 
eddy diffusion coefficients were also suggested. To achieve the 
recommended N 2 densities as well as to incorporate some recommended 
revisions in Reber' s temperature-height profile, in the 100-115 km 
region, it became necessary to lower the temperature in the 86 to 90 km 
region. This was accomplished by a revision of the model of Task 
Group I through coordination with that group. 

At the December meeting the n th version of the Reber model was 
reviewed, and further modifications were suggested. One member who had 
made no overall proposal still had misgivings, but the committee agreed 
to adopt, for submission to COESA, the Reber model of that date with 
yet a few modifications. 

Various members of the committee agreed to participate in the 
preparation of the report as follows: 

Principal Discussion of the Model -- Reber 
Composition Measurements -- Nier and Moe 

Eddy Diffusion Considerations -- Zimmerman 
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Background, and Lower Boundary Conditions 
Including discussion of Molecular Scale 
Temperature and Geopotential Minzner 

II. LOWER BOUNDARY CONDITION 

Because of the necessity for continuity with the "U.S. Standard 
Atmosphere, 1962" at a height of 50 kilometers, it was necessary for 
the several task groups to generate their respective models sequentially. 
Task Group II would therefore use the temperature and density of the 
model of Task Group I, at some point in the overlap region between 80 
and 100 km, as the lower boundary values for its model. 

The purpose of the regions of overlap for the three task groups 
was to allow needed flexibility in the generation and matching of the 
successive models. For example, the temperatures in the recommended 
isothermal region between 85 and 90 geopotential kilometers as orig- 
inally submitted by Task Group I were sufficiently high so that no rea- 
sonable temperature-height profile above 100 km would yield acceptable 
densities at 150 km. In the ensuing coordination between Task Group II 
and Task Group I it was finally agreed that the temperature gradient of 
-1.8 Kelvin degrees Der geopotential kilometer (K/km*) between 71 and 
85 km' would be changed to -2.0 K/km', thereby producing a reduction of 
2.8 K (from 189.45 K to 186.65 K) at 85 geopotential kilometers (km'), 
and simultaneously allowing for an acceptable value of mass density at 
150 km geometric kilometers (km) with a reasonable temperature profile 
between 85 km' and 150 km. (The above revision was preferred by Task 
Group I over a proposed extension of the -1.8 K/km> gradient from 85 km' 
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to 86.5 km' whereby the temperature of the isothermal layer would have 
been reduced from 189.45 to 186, 75 "K, while the base of this layer 
would have been raised 1.5 km'. 

The transition region between the temperature profile of the model 
of Task Group I and that of Task Group II also serves as a region of 
transition between two different scales of temperature, (i.e. molecular 
scale temperature T M , and kinetic temperature T) as well as between two 
different scales of height measurement (i.e., standard geopotential 
kilometers km', and geometric kilometers km). In the "U.S. Standard 
Atmosphere, 1962", the model was divided at 90 geometric kilometers 
with that part from 0 to 90 km defined in terms of T f/J and H, and the 
part above 90 km defined in terms of T and Z. A similar policy is 
being followed in the proposed revisions to the 1962 Standard. It is 
the recommendation of Task Group II, however, that the height at which 
this transition occurs be set at 86 geometric kilometers, or 84.8520 
km' where the molecular scale temperature according to the revised model 
of Task Group I is 186.9460 K. It is desirable that the transition 
occur at an integer value of geometric height measure for ease in inte- 
gration at higher altitudes in intervals of one geometric kilometer. 

Because the model being proposed by Task Group II implies the 
existance of some atomic oxygen at 86 kro, the mean molecular weight M 
applying to that height has been taken to be 28.940. At this height 
the ratio M/M q , therefore,, has the value 

(28.940/28.9644) = .9991645 

and the kinetic temperature has the value of 186.789805 which when 
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rounded to 5 significant figures is 186. 79K. Using a Task Group I value 

of 6.421062 X 10 -5 kg/m 3 for the density at 71 km where the value of 

T is 214.65 K, and assuming a constant gradient of -2K/km from 
M 

71 km to 84.8520 km or 86.000 km where the value of must therefore 

be 186.9460 K, the value of density has been computed to be 6.957839 x 

10 -6 kg/m^. The values representing the boundary conditions for the 

lower end of that portion of the Proposed Standard Atmosphere above 

86 km as recommended by Task Group II, are therefore as follows: 

H = 84.8520 km 

Z = 86 km 

M/M = .9991645 
o 

T = 186.9460 K 
M 

T = 186.789 8 K 

A brief discussion of history of the development and use of 
molecular scale temperature and geopotential is given in the Appendix. 

III. PHILOSOPHY AND CONSTRAINTS OF THE MODEL 

In addition to the basic guidelines already mentioned, it was 
felt that the form of the mathematical expressions defining the temp-’ 
erature and number- density profiles as a function of altitude should 
make the model useful as a theoretical tool throughout the altitude 
range. A number of implications follow from this philosophy: 

a. The temperature should be express able as a smooth mathematical 
function of geometric altitude, with a smooth second derivative. In 
particular, it was felt desirable to use the exponential (Bates, 1959) 
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profile in the middle and upper thermosphere as this form is well 
known, widely used, and permits the utilization of the Walker (1965; 
Bates, 1959) technique for analytically representing upper atmosphere 
number densities. 

b. The functions representing the temperature profile should be 
readily adjustable to allow approximation of varied data sets. 

c. The functions relating number densities to altitude should be 
physically meaningful and expressable analytically. 

The temperatures and gas -specie number densities of the model have 
to be consistent with inputs from a variety of sources: 

a. At the lower altitude boundary the temperature and number 
densities must match with the model recommendations of Task Group I, 

a match which is somewhat complicated by the fact that the model of TGI 
is defined in terms of geopotential altitude and molecular scale temp- 
erature, while TG II uses geometric altitude and kinetic temperature 
(see Section II). 

b. In the region between the lower boundary and about 130 km, the 
temperature and mass density profiles should match the data available, 
which come largely from rocket-borne pitot measurements, falling-sphere 
measurements and Thomson incoherent scatter measurements. However, the 
average value of the N 2 density above 150 km is fairly well established, 
and this has a large influence on the choice of temperature profiles in 
the region below this altitude. See the discussion in IV.C.l. 

c. At 150 km the composition should match the TG II recommendations 
shown in Table I and discussed in Section V. 
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d. The largest body of data available on the neutral composition 
of the thermosphere (as opposed to the larger data set available on 
total density) was that obtained from the quadrupole mass spectrometer 
on the 0G0-6 satellite (e.g. Hedin, et al., 1972). Since these data 
are referenced primarily to an altitude of 450 km it was decided to 
extend the calculations for the Task Group II model to at least that 
altitude to allow inclusion of this large and unique data set. (Table I) 

e. The decision had been made by the COESA at the September 1971 
meeting that the average conditions to be modeled could be well approxi- 
mated by using an exospheric temperature of 1000° K. 

f. At altitudes above about 150 km, the total density and its 
scale height should be consistent with the large body of data determined 
from satellite drag. 

It has been borne in mind throughout that many of the parameters 
and profiles used and calculated are dynamic by nature and any steady- 
state description is only an approximation to the true state-of-affairs . 
Examples are the temperature profile which generally exhibits wave-like 
structure, the atomic oxygen profile which calculations show to be 
extremely time dependent with a significant diurnal component, and the 
helium profile with an annual component. For the purposes of this 
model, choices were made for the best average value or profile to be 
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used or matched. 



IV. MODEL 


A. Temperature 

The temperature profile is described in four altitude ranges: 

1. z q (86 km) to z 1 (91 km): isothermal at 186.79°K. 

2. z (91 km) to z 9 (110 km): A portion of an ellipse is 
used here, assuring a smoothly monotonically increasing temperature 
with altitude, with sufficient generality to match the temperature and 


its gradient at the end points of the region of definition. 

„ , ,2 1/2 

T(z) = T + b 2 (1-^^ ) 


where T 


2 2 
YT - T + T 
1 1 2 2 1 

Y + 2T 1 -2T 2 


dt 

Y = dz 


110 (Z 2* Z 1^ ’ 


( 1 ) 


b (z 2 - z x ) 
Cb 2 -(T 2 -T c ) 2 ) 1/2 


b = CT 1 - T c ) . 


3. z 2 (110 km) to z^ (125 km) : 


dT(z) 

dz 


12°/km 

= 240 u K, and 
= 420° K. 


3 



( 2 ) 


4. Zg (125 km) to 500 km: 


TU) = T oo - ( T co - T 3 ) eX P C-° $), 
where = exospheric temperature 

a = 41 = 0.02069 km' 1 , 

T -T 7 dz ‘ 

co 3 z 3 

(Z-Z 3 ) (R + z,) 

C(z) = — , and 

R + z 


R = earth radius = 6356.77 km. 


B. Densities 


The steady-state vertical distribution of a minor gas specie of 
number density n.. and mass nr is governed by the vertical component of 

the momentum equation for that gas (e.g. Colegrove, et al., 1965): 

dn . ( 1+ a . ) n dT n . 

n.v. + D. (-r-^ + -g 1 - ) 

li l Mz T u J 


dz H. 

l 


dn. 1 dT n. 

+ K (^ +f d7 +g ±-) = 0; 


where = the flow velocity of i specie, 

th 

= local molecular diffusion coefficient for i specie 
diffusing through the major, background gas, 

= thermal diffusion factor for i 

(= -0.4 for He; = 0 for all others), 
kT 

Hi = jj— - = local scale height for i, 
i® 


k = Boltzmann constant, 
g = local acceleration of gravity. 
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K = eddy diffusion coefficient, 
kT 

H = — = local scale height of the atmosphere, and 
Mg 

M = local mean molecular mass of the atmosphere. 

1. Molecular nitrogen 

Equation 3 is used to calculate the distribution of each 
gas except molecular nitrogen; since N 2 is the major gas in the lower 
part of the thermosphere a different scheme is used to describe its 
distribution. From the lower boundary altitude of 86 km up to about 
100 km the atmosphere is well mixed, as the eddy processes dominate the 
molecular diffusion, while above about 100 km molecular diffusion 
dominates. Also, because the molecular mass of N 2 (28.0134 amu) is 
quite close to the mean mass of the atmosphere before dissociation and 
diffusion become important (28.9644) , the transition from mixing to a 
diffusive distribution has little effect on the vertical distribution 
of molecular nitrogen. This allows the use of a simplified version of 
Equation 3 for N 2> with the flow and the eddy terms set equal to 0: 

dn (N ) 1 dT n(N ) 

~ dz + T dT + HT = °' 

2 

Thus, the distribution of molecular nitrogen is affected only by the 
temperature profile. To account for the slight change from mixing to 
diffusion, the mean mass is used in the calculation from 86 km to 
100 km, and the N 2 mass is used above 100 km. 
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2. Atomic and Molecular Oxygen, Helium and Argon 
Ideally, Equation 3 is solved in conjunction with the equation 
of continuity (Colegrove, et al. , 1965; Keneshea and Zimmerman, 1970), 


t— (n. v. ) = S (z) 
dz v l l 


where S denotes production and loss terms. However, for the purpose of 
the single-profile, steady-state model being generated here, it was felt 
that a sophisticated and detailed calculation of this nature was not 
appropriate. Instead the flux terms, n^v^, is artificially adjusted to 
include the effect of photo chemical production and loss on the vertical 
distributions of atomic and molecular oxygen. For helium and argon the 
flux term represents only the vertical flow. 


Equation 3 is integrated directly to obtain 

, r z 

n. (z )T(z ) 

. . i o o r 

n (z) = exp [ - j 

T (z) X 


v. 

l 


(f(z) + D.+K ) dz], 


(4) 


where 


D. 

l 


V 1+a)+K dT 


K 


£( ^ Z - ) H. (Di + K) + T(D i +K) dz + H ^ + K) 

The eddy diffusion coefficient profile used in the calculation is of 

2 


the foim (Appendix II). 

J 


K(z) » 3.0 x 10' exp (z-88 )) j z ^ 88 km 


= 3.0 x 10 7 exp ) 2 , z > 88 km. 

The molecular diffusion coefficients represent diffusion through 
molecular nitrogen, and have the general form 
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D. (zj = “i_ 


C 


T(z) 

273.16 


) 


B. 

1 


cm 


N(z) ' ' ' sec 

where N(z) is the total number density. The coefficients a and (3 axe 

listed in Table II. The calculation of the total number density used 

for D. (z) is done sequentially, with N„ providing the N(z) for the 

calculation of the atomic and molecular oxygen diffusion coefficients. 


and n(N 0 ) + n(0)+n(0.J providing N(z) for the argon and helium coefficients. 

The flux term, ^i_ , in Equation 4 is represented by the integrable 

FT+K 

i 

expression 


= AjCz-a^ 2 exp (-b^z-a^ 3 ) + A 2 (a 2 -z) 2 exp (-b^-z) 3 ) . 
i 

The constants A , A 2 , a 1 , a 2> b^ and b 2 are determined such that appro- 
priate densities are determined at 450 km for 0 and He, and at 150 km 
for 0, 0 2 , He and Ar. The constant A 2 = 0 for all species except 
atomic oxygen; the extra term for 0 is needed to generate a maximum in 
the density profile (chosen to be at 97 km) , reflecting the increased 
loss by recombination at lower altitudes. The flux terms for 0 and 0 2 
are based on, and lead (qualitatively) to the same results as those 
derived from the much more detailed calculations by Colegrove, et al. (1965) 
and Keneshea and Zimmerman (1970). Table II lists these coefficients. 

C. Discussion 

1. N 2 density and the temperature profile 

As noted in paragraph B.l, the N 2 density at any altitude is 
sensitive primarily to the temperature profile at lower altitudes. 

This fact has serious' implications when there are a number of data sets 
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to be matched,, as in the case here, where the lower boundary conditions 
are given (at 86 km) , a density for N 2 at 150 km is given, and there 
are some temperature data available in the region between these two 
altitudes . 

The temperature data consist mainly of recent pitot -tube measurements 
(J. Theon and J. Horvath, private communication) and incoherent scatter 
measurement (e.g. R. Wand, private communication). These two data sets 
are quite consistent in one particular feature: the temperature profile 

between about 105 km and 125 km appears to have a constant gradient of 
approximately 18°/km. The proposed model does, in fact, exhibit a con- 
stant gradient in this region, but it is 12°/km, only 2/3 of the measured 
value. Attempts to incorporate higher gradients lead to unacceptably 
high values for N 2 densities above 150 km, and herein lies the dilemma. 

The lower boundary parameters and the N 2 density at 150 km reflect 
the results of many measurements in which there is a high degree of 
confidence so it is not likely that these data are in error. It is not 
clear whether the recent measurements of the temperature profile over 
estimate the gradient, or whether the three inputs are inconsistent in 
that they are not true averages over the same sets of conditions. 

2. Dynamic characteristics 

As noted earlier, most of the properties being modeled are 
time-dependent by nature, and any steady-state description has to be 
used advisedly. Examples are the diurnal photo-chemical variations in 
atomic and molecular oxygen densities and the longer tern, dynamically 
induced variations in helium and argon densities. The proposed model 
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COMPOSITION AND TOTAL DENSITY AT 120 km, 150 km and 450 km 
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includes the provision to represent deviations from diffusive equilibrium 
profiles in the middle thermosphere which are becoming more and more 
accepted as physically real. 

V. COMPOSITION 

Atmospheric densities in the 100-200 km altitude range, computed 
from composition measurements made with rocket-borne mass spectrometers, 
have always been lower than values inferred by downward extrapolation of 
drag measurements on satellites having higher altitude orbits. While 
it has been recognized that there might be some error in the drag 
coefficient upon which the drag measurements depended, the general feel- 
ing has been that the composition measurements were in error. In par- 
ticular, because of the highly reactive nature of atomic oxygen it has 
been assumed that this constituent was largely lost in mass spectrometer 
ion sources and hence grossly underestimated. Early mass spectrometric 
values such as those of Meadows and Townsend [1960] or Pokunkov [1960] 
were extremely low, undoubtedly owing to the loss of atomic oxygen on 
the extensive surfaces of their instruments. With the advent of "open" 
source instruments such as those of Schaefer [1963] and of Nier et al. 
[1964] much higher values were obtained. Even so, it was recognized that 
the losses might still be considerable. 

Hall et al [1965, 1967], using EUV extinction measurements made 
with rocket-borne UV spectrometers, found atomic oxygen abundances in the 
150-200 km altitude range to be considerably above those reported from 
rocket-borne mass spectrometer measurements. Results extrapolated 
downward from 0S0-III measurements [Hinteregger and Hall, 1969] gave 


similar results. The absolute numbers given are in some doubt, however, 
in view of the uncertainty in the absorption cross section employed 
for atomic oxygen [Moe, 1970]. 

vonZahn [1970] summarized the situation at 150 km as of 1970 making 
use of the relevant published composition values, mass spectroscopic as 
well as UV extinctions, available at the time. After evaluating the 
data he concluded that the most consistent agreement between drag and 
mass spectroscopically determined mass densities was obtained if one 
assumed that drag determined densities were high by 10% and mass spec- 
troscopically found atomic oxygen values were low by an appreciable fac- 
tor, perhaps as much as 4. Accordingly he recommended particle densi- 
ties at 150 km as follows: 

n(N 2 ) = 2.6 x 10 10 cm" 3 
nfO^] = 2.5 x 10^ cm 3 
n(Ar) = 5 x 10 7 

n(0) = 2.3 x 10 10 

-12 -3 

and mass density p = 1.96 x 10 g cm His drastic increase in the 
amount of 0 seemed justified in part by measurements made at 120 km 
with a helium-cooled rocket-borne mass spectrometer which gave an 
appreciably higher value of O/O2 at 120 km than had ever been reported 
in the literature [Offermann and von Zahn, 1971]. 

Nier [1972] ,on the other hand, pointed out that since atomic 
oxygen is a major constituent of the atmosphere in the neighborhood of 
150 km any arbitrary increase, such as by a factor of substantially more 
than two in its measured abundance relative to other constituents, would 
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destroy the excellent agreement between mass density scale heights com- 
puted from mass spectrometer composition measurements and those found from 
drag measurements on low altitude satellites such as 0V1-15 [Champion 
et al , 1970a] and 0V1-16 [Champion et al, 1970b]. He subsequently re- 
inforced his argument through laboratory experiments [Nier et al, 1972] 
and [Lake and Nier, 1973] in which it was shown that it was not likely 
that atomic oxygen densities measured with instruments such as he and 
colleagues previously used in rocket flights were low by a factor of 
more than two. 

More recently Taeusch and Carignan [197] in an extrapolation of 
0G0-6 composition and drag-determined densities down to 150 km concluded 
that the 150 km atomic oxygen value given by von Zahn [1970] and employed 
by Jacchia in his 1971 model was too high. They prefer a number about 
20 percent lower but still considerably above the average value found 
with rocket-borne mass spectrometers. Their n(N 2 ) and values at 

150 km, on the other hand, are about 25 percent higher than values gen- 
erally found with rocket-borne mass spectrometers. 

Moe [1973] completed a comprehensive study of drag measurements 
with satellites as well as of published values of atmospheric composition 
by all methods, correcting drag measurements for effects due to accommo- 
dation coefficients, and composition measurements for possible errors in 
instruments due to surface effects. His nCN^) , nCC^) and n(Ar) values 
at 150 km agree closely with those given by von Zahn [1970] , which are 
essentially the abundances found by rocket-borne mass spectrometers. 

His n(0), however, is about 20 percent lower than von Zahn's, in agreement 
with Taeusch and Carignan [1972]. 


The concentration of helium in the lower thermosphere at mid- 
latitudes is known to vary by a factor of as much as 10 between summer 
and winter. Also below 150 km it appears not to be in diffusive equil- 
ibrium. The values presented in the present report fall between the 
extremes observed in observations. 

While some of the variations reported in n(N 2 ), n(0 2 ) and n(Ar) 
measurements in the 100-200 km range are almost certainly due to errors 
in measurements, some must be attributed to true atmospheric variations. 
The values of n(N 2 ) , n(C> 2 ) , n(Ar) at 150 km used in constructing the 
present model are nominal values and are the best estimates at the pre- 
sent time. Each is believed to be correct to 25 percent. Because of 
the uncertainty in the amount of atomic oxygen lost in rocket-borne mass 
spectrometers, the value of n(0) at 150 km is based in part on values 
extrapolated downward from measurements made at higher altitudes with 
satellite-borne instruments in which it is believed that the atomic oxygen 
loss can be properly evaluated [Hedin et al, 1973] and in part on mass 
densities found from satellite drag and corrected for the other con- 
stituents (N 2 , 0 2 and Ar) , which can be measured accurately. It appears 
likely that the n(0) values given in the table are maximum values, as 
they are based on the assumption that atomic oxygen is strongly absorbed 
in mass spectrometers used in rocket studies of the lower thermosphere. 
This view may be too pessimistic but it does not seem probable that 
values given could be high by a factor as large as two. 
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APPENDIX I 


MOLECULAR SCALE TEMPERATURE AND GEOPOTENTIAL 

The concept of a "derived 'Temperature' based upon an assumed con- 
stant mean molecular weight y of 28.966 g/mole for the atmosphere" was 
first applied to atmospheric models in a paper by the Rocket Panel 
(1952) .. 

The lack of a specific name for this derived "temperature led to 
the adoption of the name molecular scale temperature in the ARDC Model 
Atmosphere, 1956 (Minzner, 1956). This concept was carried into the U.S. 
Extension to the ICAO Standard Atmosphere (Minzner et al., 1958), and 
into the U.S. Standard Atmosphere, 1962. 

The concept of geopotential as a. measure of height stems from 
V. Bjerknes et al., (1910), who "made use of the term 'dynamic height' in 
referring to the geopotential of a point because the latter is preferable 
to geometric height in meteorology as a representation of the vertical 
coordinate of the point" and who "proposed as the unit of geopotential the 
so-called geodynamic meter (gdm.) or dynamic meter for short" (list, R. J., 
1951). With a slight revision in definition, this concept led to the 
geopotential meter (List, R. J., 1951) and to the standard geopotential 
meter, the latter of which was first applied explicitly to standard 
atmospheres in NACA Report 1235, (Annon. , 1955). The use of this con- 
cept was continued in the ARDC Model Atmosphere, 1956 (Minzner, 1956) , in 
the U. S. Extension to the ICAO Standard Atmosphere (Minzner, et al., 

1958) as well as in the U.S. Standard Atmosphere, 1962. The standard 

34 



geopotential meter is that unit of geopotential defined in terms of the 

_2 

standard seal- level value of the acceleration of gravity, 9,80665 m sec . 
Implicitly, geopotential has been used as a measure of height in all 
earlier U.S. Standard Atmospheres in which the tabulated values of 
pressure were calculated on the basis of a value of the acceleration of 
gravity, which is invariant with height. 

The concept of geopotential has some very meaningful uses in 
meteorology. In Standard Atmospheres, however, the use of this concept 
was introduced, either implicitly or explicitly, prior to the development 
of high-speed digital computers, as a means of avoiding the more compli- 
cated equations which result when the pressure-height relationship is 
developed in terms of an acceleration of gravity varying as some function 
of height. The use of geopotential has been maintained in all standard 
atmospheres to date in order to avoid any revision of the lower portion 
of the tables which have represented the established standard for the 
past 20 to 50 years. 

The introduction of the concept of molecular scale temperature T M 
came at a time when standard- atmosphere tables were being extended to 
heights where the composition, and hence the mean molecular weight were 
unknown. The use of T M not only avoided the problem of determining or 
guessing at a value of mean molecular weight at high -.altitudes, but 
also eased the problem of hand or desk-computer calculations by leading 
to simpler equations than would have resulted if specific height func- 
tions had been introduced for both kinetic temperature and mean molecular 
weight. 35 


In an attempt to further simplify the calculation of standard 
atmospheres before th© wide use of these high-speed computers, Brombacker 
(1953) attempted to combine the height- dependent acceleration of gravity 
the height-dependent varying molecular weight and the height-dependent 
kinetic temperature into a single variable which he called scale-height 
temperature. This concept, however, was never adopted in standard 
atmospheres . 

In each of the following; the ARDC Model Atmosphere, 1956; the 

U.S. Extension to the I.C.A.O. Standard Atmosphere; and the ARDC Model 

Atmosphere 1959 (Minzner, 1959); the entire model was defined in terms 

of molecular-scale temperature in units of Kelvin degrees K, and 

geopotential H in units of km . Since most aeronemists were unfamiliar 

with these quantities and preferred kinetic temperature T, and geometric 

heights Z it was decided that in the U.S. Standard Atmosphere 1962, the 

model was divided at 90 geometric kilometers with that part from 0 to 90 

km defined in terms of T.. and H, and the part above 90 km defined in 

M 

terms of T and Z. 
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APPENDIX II 


COMPLETE TIME DEPENDENT CALCULATION 

The philosophy of these calculations is tne generation of self consistent 
Nitrogen, Oxygen, Ozone and Argon diurnally averaged atmospheres that fall 
within the species concentration constraints and altitudes recommended oy 
Task Group II. The species profiles are the result of a time dependent 
photochemical calculation using molecular and turbulent transport which 
incorporates the latest measured chemical reaction rates, temperature, 
solar flux and turbulent diffusion coefficients into the coupled equations of 
motion and continuity. 

The species considered are O^, O O^, O^D, O 2 -4 g. OH, H, HO 2 . 

H O, H O H-> and Ar. The diurnally varying number densities, from 50 
2 2 2 ’ ^ 

to 150 km, are obtained through a finite difference solution of a system of 
mass and momentum conservation equations (Shimazaki, 1967, Keneshea ana. Z 
merman, 1970). The numerical approach is essentially that introduced 
by Shimazaki but modified at the boundaries and in the volume integrations, 
following George et al (1972). The chemical production and loss rates for 
each species are displayed in the reaction scheme of table 1. The intensity 
of the solar flux is that reported in Ackerman (1 970)(figure 1) and the 
absorption cross sections are those of Hudson (1972). The temperature 
profile and the mean molecular mass used up to the turbopause are those 
recommended by Task Group II. Using these data, the initial species 
distributions are calculated assuming mixing to the turbopause and 


diffusive equilibrium above. 
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The total number density is obtained by integrating the hydrostatic 
equation where the sea level mass density and the mean molecular weight 
profiles are taken from the U.S. Standard Atmosphere (1962). 

Beginning with the static profiles described above, a steady state 
solution is determined, which in effect significantly reduces computer 
time to the final output. The solutions are then continued for ~ 10 days 
using a fully implicit finite difference technique, a variable time step, 
limited by species changes of 2 efold and a height step of 100 m. This 
stringent height step was shown to be necessary to restrict the errors 
generated by species gradients when height steps larger than 100 m 
were used. 

The turbulent diffusion coefficients used are based upon observations 
of turbulence in chemical trails, the observed altitudes of cessation and 
the measured turbulent diffusion coefficients (Zimmerman et al, 1970; 
Keneshea and Zimmerman, 1970). The turbopause selected was 
strongly influenced by mass spectrometer measurements (Van Zahn, 1970) 
where, in general, the reported altitude of transition from mixiiig to 
diffusive equilibrium was ~ 102 km. Thus the description of the turbulent 
diffusion coefficients is given by 

2 

K = Aj exp (- 2. 303 
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where the values selected are 



Z = 
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88 km, 
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a - < 

1 25 km, Z 

£ 88 

and 

A = 

3 x 10 7 . 



The lowest altitude ( ~ 50 km) value determined by the above relation 
approximately matches the reported measurement (K = 10 cm /sec) 

by Beaudoin et al (1967). The peak value used falls within the 
experimentally determined values of the vertical turbulent diffusion 
coefficients (Philbrick et al. , 1973), (Keneshea and Zimmerman, 1970). 
Because of the lack of chemical tracer wind and turbulence measurements 
in the altitude region 60 to 90 km, we are forced to assume an exponential 
fit between the measurements at 50 km and 90 km. 

RESULTS 

The time dependent calculations are continued for ~ 10 days until 
the species concentrations reproduce themselves within 1% over a 
diurnal cycle, or what is phrased as arriving at diurnal reproducability. 
The diurnal average of [O], and is then calculated and 

extrapolated to 250 km by assuming diffusive equilibrium above the upper 
boundary. Figures 1 and 2 show the initial conditions of [N 2 ], 
Temperature, and the eddy diffusion coefficient used in these one 
dimensional calculations. Figure 3 shows the diurnally averaged species 
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profiles [O], [0 2 ], [O^], [H] and [Ar], with the recommended values 
of Task Group II at 150 km, also shown. For the turbopause value of 
~ 102 km (defined here as the altitude at which the molecular self 
diffusion coefficient is equal to the turbulent diffusion coefficient) the 
theoretical values show quite good agreement with the above recommended 
values at 150 km. Figure 4 presents the argon to nitrogen and the O to C> 2 
ratios. Observe the gradual separation of argon from the mixed region 
below the turbopause to completely diffusive separation some, distance 
above the turbopause. Because of this fairly large transition region, a 
simple extrapolation downward of an argon measurement until the ground 
layer mixing ratio is achieved, can be misleading when trying to describe 
the turbopause height. This effect, while not shown here, is even 
greater for helium. 

Thus in conclusion, it has been demonstrated that an internally 
self consistent picture of the density structure of the upper mesosphere 
and lower thermosphere may be achieved using measured values of 
solar flux, reaction rates and/or derivatives of measured vertical 
turbulent transport parameters deduced from chemical trail studies. 

The fairly good agreement with the limits of these species at 150 km 
placed upon them by Task Group II may be fortuitous, particularly in 
the light of the gross differences and uncertainties in the oxygen mass 
spectrometer measurements. However, it does show that even with 
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these uncertainties the physical insight and intuition of Task Group II, 

in estimating the above quantities is exceedingly good. 
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